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It is by now well established that the analysis of the 
level statistics is one of the main tools in the study of 
quantum complex systems. Moreover the spectrum, un- 
like the wavef unctions, is easily accessible either numeri- 
cally or experimentally. Part of this interest stems from 
the fact that, once the model-dependent spectral den- 
sity is extracted from the spectrum, the level correla- 
tions of apparently unrelated models shows striking uni- 
versal features in a variety of physical situations. For 
instance, the celebrated Bohigas-Giannoni-Schmit con- 
jecture (BGS) states the level statistics of a quantum 
system whose classical counterpart is deterministic but 
fully chaotic does not depend on the microscopic details 
of the Hamiltonian but only on the global symmetries of 
the system and coincides with those of a random matrix 
with the same symmetry (usually referred to as Wigner- 
Dyson statistics (WD)0]). 

Remarkably the same WD also describes Q the spec- 
tral correlations of a disordered system Q in the metallic 
limit. By contrast, for disorder strong enough, eigen- 
states localization becomes important, the spectrum is 
not correlated and the level statistics is universally de- 
scribed by Poisson statistics. For deterministic systems 
Poisson statistics appears provided that the classical dy- 
namics is integrable [sj. 

Despite its robustness, the universality associated to 
WD has also clear limits of applicability. For instance, 
the quantum properties of Hamiltonians whose classical 
phase space is a superposition of chaotic and integrable 
parts are supposed to depend dramatically on the details 
of the Hamiltonian and consequently their properties are 
non-universal. Similarly, for finite disordered systems 
in the metallic regime, the dimensionless conductance 
g — Ec/ A {Ec, the Thouless energy, is a scale of en- 
ergy associated with the classical diffusion time through 
sample and A is the mean level spacing) sets the number 



of eigenvalues whose spectral correlations are universally 
described by WD. 

Universality in the spectral correlations has also a 
counterpart in the eigenfunctions properties. Thus Pois- 
son statistics is associated with exponential localization 
of the eigenfunctions and WD is typical of systems in 
which the eigenstates are delocalized through the sample 
and can be effectively represented by a superposition of 
plane waves with random phases. 

Recently it was reported H, 0, |1| that the quan- 
tum properties of a disordered conductor at the metal- 
insulator transition (usually referred to as the Anderson 
transition (AT)) are to certain extent also universal. Fea- 
tures related to this new universality class include multi- 
fractal eigenstates and level statistics given by critical 
statistics '7](see below for a definition). Intuitively 
multifractality means that eigenstates have structures at 
all scales. Roughly speaking the amplitude of probabil- 
ity of a multifractal eigenstate has peaks ('probability 
splashes') at all scales decaying as a power-law from its 
maximum. Consider the volume of the subset of a box 
for which the absolute value of the wave function is 
larger than a fixed number M. If this volume scales as 
L'^ (with d* < d), then d* is called the fractal dimension 
of In case the fractal dimension depends on the value 
of M , the wave function is said to be multifractal. On a 
more formal level multifractality is defined either through 
the box counting method (see [2^ and the eigenvector 
analysis section below) or the anomalous scaling of the 
eigenfunction moments Vq = J d''r|'i/'(r)P'' with respect 
to the sample size L as Vq oc L~^i^'^~^\ where Dq is a 
set of exponents describing the AT 0. 'Critical statis- 
tics' 0,13 (tti^ level statistics at the AT) is intermediate 
between Wigner-Dyson and Poisson statistics. Typical 
features include: scale invariant spectrum Q, level re- 
pulsion and asymptotically linear number variance Q. 
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Level and eigenfunction correlations at the AT are said 
to be universal in the sense that typical features of criti- 
cal statistics as level repulsion, scale invariance and linear 
number variance (E^(L) = (L^) — (L)^ = for x < 1 
and L ^ 1) do not depend on boundary conditions, the 
shape of the system or the microscopic details of the dis- 
ordered potential Tl'l . However the slope of the number 
variance or the functional form of certain level correlators 
as the level spacing distribution P{s) (the probability of 
having two eigenvalues at a distance s) may depend on 
additional parameters as the dimensionality of the space. 

Another argument reinforcing the universality of criti- 
cal statistics is the fact that, as for WD, which describes 
the level statistics of a Gaussian random matrix model, 
critical statistics has also been found in a variety of gen- 
eralized random matrix models: based on soft confi ning 
potentials , effective eigenvalue distributions 0, Q| 
related to the Calogero-Sutherland model at finite tem- 
perature and random banded matrices with power-law 
decay . The latter is specially interesting since an AT 
has been analytically established by mapping the prob- 
lem onto a non linear a model. 

Finally we recall that critical statistics is not related 
to any ergodic limit of the quantum motion as in the 
case of WD. Consequently it is capable to describes non 
trivial dynamical features and its limit of validity is not 
restricted to a scale given by the dimensionless conduc- 
tance g (as for WD). 

A natural question to ask is whether critical statisti- 
cal and multifractal wavefunctions are exclusive of disor- 
dered systems or may also appear in deterministic quan- 
tum systems. Indeed, in a recent letter [T^ we have es- 
tablished a novel relation between the presence of anoma- 
lous diffusion in the classical dynamics, the singularities 
of a classically chaotic potential and the power-law lo- 
calization of the quantum eigenstates. Specifically, for 
a certain kind of singularity (log for 1+lD system) as- 
sociated to classical 1// noise, it is found that the level 
statistics is described by critical statistics and the eigen- 
states are multifractal with a Dq quantitatively similar 
to the one at the AT. These results are universal in the 
sense that neither the classical nor the quantum proper- 
ties depend on the details of the potential but only on 
the type of singularity. 

Other non-random systems whose level statistics has 
been reported to be similar to critical statistics include: 
Coulomb biUiards 0|, Anisotropic Kepler problem [Tsjl . 
generalized kicked rotors ^3 '^'^^ pseudointegrable bil- 
liards [2(1 l2lj | . For the latter the dynamics is interme- 
diate between chaotic and integrable and, in order to 
fit the spectrum, Bogomolny and coworkers introduced 
|2fl| | a purely phenomenological short range plasma model 
whose joint distribution of eigenvalues 2QJ is given by the 
classical Dyson gas with the logarithmic pairwise inter- 
action restricted to a finite number k of nearest neigh- 
bors. Explicit analytical solutions are available for gen- 
eral k. For instance for k = 2 (usually referred to as Semi- 
Poisson statistics SP), i?2(s) = 1 - e'^', P{s) = ise''^'' 



and T,'^{L) = L/2 + (1 - e-^^)/8 where i?2(s) is the two 
level correlation function (TLCF). In passing we mention 
that SP can also be obtained by removing every k eigen- 
value out of a spectrum with Poisson statistics. It turns 
out that this short-range plasma model reproduces typ- 
ical characteristics of critical statistics as level repulsion 
and linear number variance with a slope depending on k. 

However SP is quantitatively different from critical 
statistics. In critical statistics, as mentioned above, the 
joint distribution of eigenvalues can be considered as an 
ensemble of free particles at finite temperature with a 
nontrivial statistical interaction. The statistical interac- 
tion resembles the Vandermonde determinant and the ef- 
fect of a finite temperature is to suppress the correlations 
of distant eig envalues. In the case of SP this suppression 
is abrupt [20j since only nearest neighbor levels can inter- 
act. Thus critical statistics and SP share similar generic 
features but are in principle quantitatively different. 

The aim of the letter is to establish under what generic 
circumstances one may expect SP in the context of non- 
random Hamiltonians. 

For the sake of clearness we enunciate our main con- 
clusions: We have found that the appearance of SP can 
indeed be traced back to the presence of a certain kind of 
singularity (different from the one associated to critical 
statistics) in the classical potential. It is shown that SP 
is indeed robust under arbitrary smooth perturbations of 
the classical potential or the insertion of a magnetic flux 
provided that the classical singularity is preserved. The 
eigenfunctions associated to SP are found to be multifrac- 
tal but quantitatively different from those of a disordered 
conductor at the AT. Finally we argue that quantum 
anomalous diffusion induced by the classical singularity 
may be verified experimentally by using ultra-cold atoms 
techniques. 

The organization of the paper is as follows, in the next 
section we introduce the model: a generalized kicked ro- 
tor in a potential with a step-like singularity. In section 
two we discuss analytical results available for our model, 
then we investigate the level statistics and perform a mul- 
tifractal analysis of the eigenstates. Finally, in section 
three we examine the quantum diffusion of our model 
and discuss possible ways of experimental verification. 

I. THE MODEL 

We investigate a generalized kicked rotor in 1-1- ID with 
a step-like singularity, 

n = ^ + Viq)J2S{t^nT) (1) 

n 

with q G TT, tt), V{q) is an arbitrary nonanalytical func- 
tion with a step-like singularity. The simplest case cor- 
responds to, 

W„)^/^o, if<ze[-a,a) .3) 
1^ 0, otherwise. 
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where a sets the size of the step and ■^o the height. We 
shah see in a broad range of parameters our results do 
not depend on the specific form of V{q) but only on the 
presence of the singularity. 

The quantum dynamics is governed by the quantum 
evolution operator U over a period T. Thus, after a pe- 
riod T, an initial state tpo evolves to ^]{T) = U^o — 

e 4h e h e if^Q where p and q stands for the usual 
momentum and position operator. Our aim is to solve 
the eigenvalue problem U'^n = e~*''"/^\l/„ where is 
an eigenstate of U with quasi-eigenenergy k„. In or- 
der to proceed we first express the evolution operator 
{m\U\n) = Unrn in a basis of momentum eigenstates 
(In) = 4=1 with n = 0, . . . TV -> oo, 

Urnn = ^ J dQe'?(™-«)-*^(9) A (3) 

We remark that in this representation, referred to as 
'cylinder representation', the resulting matrix Unm is uni- 
tary exclusively in the — s- oo limit. For practical cal- 
culations this is certainly a disadvantage since besides 
typical finite size one has also to face truncation effects, 
namely, the integral of the density of probability is not 
exactly the unity and eigenvalues are not pure phases 
(e~*^" ) as expected in a Unitary matrix. 

An alternative procedure, referred to as 'torus repre- 
sentation', specially suited for numerical calculations, is 
to make the Hilbert space finite (i.e. m, n = 1, . . . TV) but 
still keep the matrix Unitary. This can be achieved by 
imposing periodicity (with period P N) in momentum 
space. With these choices h is set to the unity since the 
period in momentum P ~ iV and the period in position 
Q ^ 2tt are related by PQ = iirNh. Thus the period on 
momentum is large (P = N) and the phase space effec- 
tively resembles a cylinder (for more details see p^V 

In order to keep the kinetic term of the evolution ma- 
trix also periodic, we take T = 2ttM/N with M an in- 
teger (not a divisor of N). For technical reasons (see 
below) M is chosen to make T roughly constant (« 0.1 
in this paper) such that M/N be a good approximation 
to a irrational for every N used. The resulting evolution 
matrix (for N odd) then reads 

{m\U\n) = le-'w('"'+»') ^ (4) 

where (l){l,m,n) = 2TT{l+9Q){m~n)/N~V{2TT{l+ea)/N), 
I = -{N - l)/2, . . . (iV - l)/2 and < 6I0 < 1; 6I0 is a 
parameter depending on the boundary conditions (^o = 
for periodic boundary conditions). 

In this paper we shall use both representations depend- 
ing on the issue to be discussed. Thus for the analytical 
analysis the cylinder representation is more appropriate: 
the limit N ^ 00 can be effectively taken and conse- 
quently truncation effects are absent. For the numerical 
calculations we use the torus representation due to the 
difficulty to deal with truncation effects. An exception 



is the case of quantum diffusion where these effects can 
be accurately detected and subtracted from the calcula- 
tion. As a final remark we mention that the numerical 
evaluation of the eigenvalues and eigenvectors of U (ei- 
ther for the torus or the cylinder) is carried out by using 
standard diagonalization techniques for volumes ranging 
from iV = 500 to TV = 8000. For 6*0 = 0, parity is a good 
quantum number and consequently states with different 
parity must be treated separately. 

II. RESULTS 

For the sake of clearness we first summarize our main 
results: 

1. The level statistics associated to the evolution ma- 
trix of the Hamiltonian Eq. is scale invariant and given 
by SP in the region S = tan(uo/2) ^ 1 where vq is the 
height of the step-like potential. These results are uni- 
versal: they neither depend on the specific form of V{q) 
(provided that the singularity is preserved) nor on any 
source (as a magnetic flux) of time reversal symmetry. 

2. In the limit S = tan(wo/2) <^ 1 the spectrum is 
scale invariant and well described by Poisson statistics. 
For intermediate values of S we observe a transition from 
SP to Poisson as S is decreased. 

3. In the SP region the eigenstates are multifractal but 
with a multifractal spectrum Dq ^ A/q + Doo essentially 
different from the one observed at the AT. 

4. The classical singularity induces quantum anoma- 
lous diffusion. 

A. Analytical results 

In this section we investigate what kind of analyti- 
cal information can be obtained from the Hamiltonian 
Eq.(Q. Our aim is to show that in a certain range of 
parameter the level statistics associated to Eq. (QJ is ex- 
actly given by SP. In a first step we map, following the 
standard prescription of Ref. |23 | , the evolution matrix 
Eq.Q) onto a ID Anderson model. For smooth poten- 
tials such mapping results in a standard ID Anderson 
model with short range hopping. As known, for this kind 
of models eigenstates are exponentially localized for any 
amount of disorder and consequently dynamical localiza- 
tion has been reported |2^ in the associated kicked-rotor 
system. 

However in our case the situation is different. The clas- 
sical step-like singularity induces long-range correlations 
in the associated ID Anderson model, 

H-iAn = eni>n + ^ F{m - Tl)^^ (5) 
m 

where e„ ^ tan(rri^). 



F[m — n) = 




d6ltan(l/(6l)/2)e-^''(™-") 
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sm"/(m — n) 



= A 



A ^ tan(uo/2) and ^ — a. For e„ a random number 
from a box distribution [—W/2,W2], it has been shown 
recently (see discussion below) without invoking any 
ensemble average that the level statistics is exactly given 
by SP in the limit A 3> W. However as A gets compa- 
rable to W a shift to Poisson is observed. Thus for our 
Hamiltonian Eq.QJ the region in which SP holds corre- 
sponds with vq ^ TT {A — tan(wo/2) 1). Moreover T 
must be such that e„ ^ tan(rn^) is approximately ran- 
dom, this always occurs for T ^ 1 irrational. In the torus 
representation of the evolution matrix T must be a ra- 
tional so in order to overcome this problem we set T to 
be a good approximation to an irrational number. 
For the sake of completeness the rest of this section is 
devoted, following Ref. 24], to a technical account of the 
reasons for the appearance of SP in the Hamiltonian Eq. 
(jSJ . We recall that throughout the demonstration the di- 
agonal energy e„ is assumed to be strictly random. First 
of all we express the Hamiltonian (jSJ in Fourier space as. 



where 



and 



We fix 7 = 7r/2 (our findings do not depend on 7). 
After a simple calculation we found that Ek is not a 
smooth function (this step-like singularity is indeed the 
seed for the appearance of SP), Ek — An/ 2 for fc < tt 
and Ek = —Air/ 2 for k > n. There are thus only two 
possible values of the energy separated by a gap S — An. 
Upon adding a weak {A ^ W) disordered potential this 
degeneracy is lifted and the spectrum is composed of two 
separate bands of size ~ W around each of the bare 
points —Att/2., Att/2. Since the Hamiltonian is invariant 
under the transformation A —A, the spectrum must 
also posses that symmetry. That means that, to leading 
order in A (neglecting 1/A corrections), the number of 
independent eigenvalues of Eq. [Sjis n/2 instead of n. 

We now show how this degeneracy affects the roots 
(eigenvalues) of the characteristic polynomial P{t) — 
det{H - tl). Let 

Pdis[t) = flo + ait + . . . 

be the characteristic polynomial associated with the dis- 
ordered part of the Hamiltonian. We remark that despite 
of its complicated form, its roots, by definition, are ran- 
dom numbers with a box distribution [— W,T/F]. On the 
other hand, in the clean case. 



(tt factors are not considered). Due to the A ^ —A sym- 
metry, the full (Eq. case Pfuii corresponds with P^is 
but replacing factors by a combination {t — A)*'^{t + 
A)^'^ with ki + k2 — k. The roots of PfuU will be in gen- 
eral complicated functions of A. However in the limit of 
interest, A ^ W 00, an analytical evaluation is possi- 
ble. By setting t — ti — A we look for roots ti of order 
the unity in the A band. We next perform an expansion 
of the characteristics polynomial Pfuii to leading order 
in A. Thus we keep terms A^^^ and neglect lower powers 
in A. The resulting P/„n is given by, 

„ J- n/2-2 „ J- n/2-3 

Pfuii = h + ^ + 



n + 2 



,(i) = (t-yl)"/2(t + yl)»/2 



where the coefficients a„ are the same than those of Pdia 
above but only n/2 of them appear in the full case. The 
eigenvalues of Eq.© around the A band are = A+/3i 
with f3i a root of Pfuii ■ The effect of the long range in- 
teraction is just to remove all the terms with coefficients 
ao to a„/2 from the characteristic polynomial of the di- 
agonal disordered case. The spectrum is thus that of a 
pure diagonal disorder where half of the eigenvalues have 
been removed. The remaining eigenvalues are still sym- 
metrically distributed (the ones with largest modulus are 

well approximated by tmax — ^\J'^^f^) around A. That 
means, by symmetry considerations, that the removed 
ones must be either the odd or the even ones. This is pre- 
cisely the definition of semi-Poisson statistics. In conclu- 
sion, our model reproduces exactly the mechanism which 
is utilized in the very definition of SP. We finally mention 
that the only effect of the coefficients 3, 4 . . . , n -f 2/2 is 
to renormalize the effective size of the spectrum: ~ 2W 
for diagonal disorder and ~ 2W/^/2> for the Eq. 



B. Level statistics 

The above analytical arguments have been fully cor- 
roborated by numerical calculations. As mentioned pre- 
viously in order to avoid truncation effects we have cho- 
sen the torus representation of the evolution matrix. Our 
first goal is to show that the level statistics of our model 
with the potential of Eq. l^Jl is given by SP in the region 
S — tan(wo/2) ^ 1 for any a. In Fig. 1 we plot the 
level spacing distribution P{s) for different values of a 
and S* ^ 1. As observed, the agreement with SP is ex- 
cellent in all cases even in the tail of the distribution and 
it is not restricted to short range correlators. As shown 
in Fig. lb, the study of long range correlators as the 
spectral rigidity A3(L) = /^f (i^ - 2L'^r + r^)T?{r)dr 
further confirms this point. Up to scales of 30 — 40 
eigenvalues deviations from SP are almost indistinguish- 
able for different values of the parameters. Deviations 
for larger scales are due to well known finite size ef- 
fects. Although not shown, we have also checked that in 
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FIG. 1: (Color online) Level spacing distribution P(s) (a) and 
spectral rigidity A3 (L) (b) for the spectrum of the evolution 
matrix Eq.l^J with potential Eq.l|5J for different widths a and 
vq = TT. In all cases the size of the evolution matrix is A'^ = 
6397. The agreement with SP is excellent. 



the range of volumes accessible to nunierical techniques, 
the level statistics was to great extent scale invariant, 
namely, it was independent from the system size. Once 
the region in which SP holds has been established, we 
investigate the robustness/universality of these results. 
In order to proceed we have repeated the level statis- 
tics analysis for a potential with a non-analytical step- 
like form but perturbed by a smooth chaotic potential. 
As shown in Fig. 2, neither short nor long range spec- 
tral correlators are affected by the chaotic perturbation 
provided that the step-like non-analyticity is preserved. 
The perturbation potential Vper used in Fig. 2 is defined 
by Vper = ci cos(g) -I- C2 sin((7) + C3 cos(2q) -I- C4 sin(2(j) 
with Ci randomly chosen from a uniform distribution 
(0, 1). We have also checked that higher frequency com- 
ponents cos(3(7) do not change the results provided 

that \Vpert \ < Vq. 

In addition, we have also found that the level statistics 
is not modified if time reversal invariance is broken by 
adding a magnetic flux to Eq. (which is equivalent 
to setting T an irrational multiple of 27r; see 0]). These 
results indicate that SP is universal; namely, it does not 
depend on the details of the potential provided that the 



FIG. 2: (Color online) Level spacing distribution P{s) (a) 
and spectral rigidity A3(L) (b) of the spectrum of the evo- 
lution matrix {N = 6397) Eq.(^ with potential Eq.(|2J and 
vq = Ti and a — n/2. As observed the level statistics is not af- 
fected neither by the breaking of the time reversal invariance 
(dashed-dotted) line nor by the addition of a smooth pertur- 
bation (solid line) Vper = ci cos(g) + C2 sin(g) -I- ca cos(2g) -|- 
C4sin(2g) with random Ci (see text). 

step-like singularity is still present. 

Finally we study the transition to Poisson statistics as 
the parameter S = tan(t;o/2) goes from S ^ I (SP re- 
gion) to 5 <C 1. As shown in Fig. 3 the level statistics 
(both P{s) and A3(L)) seems to move smoothly from 
SP to Poisson. All of the typical features of criticality 
above mentioned are maintained through the transition. 
Of course parameters as the slope of the number variance 
X runs from x = 1/2 (SP) to x = 1 (Poisson) as S is de- 
creased. However, we stress that, unlike the SP region, 
it is hard to unambiguously define a universality class 
in the transition region. The point is that, in essence, 
the appearance of SP can be traced back (see analytical 
analysis) to the gap in the spectrum of the Hamiltonian. 
The step-like singularity thus separates the spectrum in 
two different bands around the only two eigenvalues in 
the clean case. As the disorder strength becomes com- 
parable with the numerical value of the bare eigenvalues 
both bands mix up and the spectrum deviates from SP. 
We could not find any sign of universality in the way 
in which this mixing occurs since it may depend on the 
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FIG. 3: (Color online) Level spacing distribution (a) and spec- 
tral rigidity (b) of the evolution matrix (TV = 6397) Eq.Q 
for the potential Eq|21for different t;o and a = 1. A transition 
from SP to Poisson is found as uo goes from tt to 0. 



microscopic details of the potential. Due to this lack of 
universality in the rest of the paper we will focus our in- 
vestigation mainly on the region of parameters associated 
to SP. 



C. Eigenvector analysis 

We now investigate the eigenvector properties of the 
Hamiltonian Eq. ^ with potential Eq. Q in the region 
5^1 corresponding to SP. 

This choice stems from the fact that, as shown above, 
the quantum properties do not depend on the details 
V{(i) but only on the presence of the step-like singularity. 
For the numerical calculation we have again utilized the 
torus representation of the evolution matrix in order to 
avoid leaking of probability due to truncation effects. 

We have two clear objectives in this section. On the 
one hand we would like to investigate whether the eigen- 
vectors associated to SP are multifractal as at the AT. 
Once this question is answered affirmatively our inten- 
tion is to provide a careful and detailed analysis of the 
anomalous dimension Dq controlling the eigenstates mul- 
tifractality. Based on the numerical findings we conjec- 
ture the relation Dq = A/q + Doc which we claim to be 
valid for all systems with SP. 
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FIG. 4: (Color online) System size dependence of the mul- 
tifractal dimension {D2(N)). For each A'^ the average is 
taken over around 36000 eigenstates of the evolution op- 
erator Eq.(|3J with the potential given by Eq. Q and 
a — 1, vo — IT. The best fitting function is a exponential 
(0.0493 ±0.0007)e-^/(22^''±i22) +D2 with D2 = 0.381 ±0.001 
(solid curve). 



To start with we give a detailed account on how Dq was 
computed. We use the standard box-counting procedure 
[25|. In doing so, for a given eigenvector -0 = X^fcLi V'fel^) 
of the system we first distribute all the components into 
Ni = N/l boxes of the same size I, then associate a prob- 
ability Pi = IV'feP to each box with k g [il, {i + 1)1); 
i = l,...,Ni. The normalized qth moments oi pi define 
the generalized fractal dimensions [26l| : 



D,iN) 



^lim^^S^ 

q - 1 s^n 



1 Pi 



log (5 



(6) 



where S = l/N is the size of the box normalized by the 
system size N . In practical calculations Dq is evaluated 
over an appropriate range of 5 by performing a linear re- 
gression of log ^^.^ with log (5 (usually a/A^ <5<l/2 
where a is a characteristic microscopic length scale of 
the system jl^)- our case a good linear depen- 
dence of logX^iJiP? log<^ is observed in the region 
5 e (0.022, 0.45) for a broad range of the potential param- 
eters , a of Eq. (0) . Dq [N) is thus computed by linear 
regression over this J- window for each eigenvector. Since 
Dq{N) is a self averaging quantity the mean value 
{Dq{N)) over all eigenstates provides with an meaningful 
description of the multifractals properties of the system. 
We observe that, even after this averaging, the statistical 
fluctuations of {Dq{N)) are still quite strong. In order 
to minimize these fluctuations we perform an additional 
average over realizations of the evolution operator Eq. 
Q until we obtain around of 36000 eigenvalues for each 
value of A^ (see Fig. 4). We generate different realiza- 
tion of the evolution matrix by varying of the period T. 
In the time reversal invariant case (T = 2ttM/N) this 
is done by picking up different values of M and in the 
broken time reversal case (T = 2t: f3 with /? irrational) by 
choosing different values of /3. 
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FIG. 5: (Color online) Generalized fractal dimensions Dq for 
the eigenvectors of the evolution matrix Eq. I^J with potential 
Eq. and a = 1, vo = TT. Solid dots are the results from Eq. 
® and the open dots are those obtained from the scaling of 
the averaged \ogVq with logA'^ (see text). The solid curve is 
the conjecture A/q + with A = 0.286 ± 0.001 and = 
0.238 ± 0.001. Inset shows the same results together with the 
prediction of critical statistics Eq. ((TJ for 6 = 0.06, 0.18 and 
0.30 respectively. 



As observed in Fig. 4, the {Dq{N)) thus obtained has 
a significant dependence on the system size which gets 
smaller as we approach the thermodynamic Umit. In or- 
der to remove this finite size efi^ect from {Dq{N)) we have 
tried different fittings, {Dq{N)) = + f{N,a,b) with 
f{N, a,b) as N ^ oo and a, 6, Dq fitting parameters. 
Thus the parameter Dq corresponds with the N oo 
limit of {Dq{N)). We found that for various values of 
q the choice of / that best describes the finite size cor- 
rections is f{N,a,b) — ae~^/^ (see Fig. 4 for the case 
9 = 2). 

After the technical analysis we are ready to present our 
results. First of all we have found that the eigenstates 



are indeed multifractal. As shown in Fig. 5, Dq (the 
— > oo of {Dq{N))) depends clearly on q. We have also 
found that Dq is very robust. Thus it is not affected by 
different choices of a, vq ^ tt, smooth perturbations to 
the potential Eq. ||2Jl, or by the breaking of the original 
time reversal symmetry. 

The next task we face is try to find an explicit expres- 
sion for Dq. Unfortunately we are not aware of any an- 
alytical technique capable of obtaining Dq exactly. The 
comparatively short range of q values for which numeri- 
cal calculations are available together with the statistical 
fiuctuations make a numerical fitting ambiguous in the 
sense that many different fitting may yield quite satisfac- 
tory results. On the other hand we think it is extremely 
important to have a proper characterization of Dq in or- 
der to completely define the universality class associated 
to SP. We found that the simplest expression for Dq still 
compatible with the numerical analysis is Dq = A/q+Dao 
(with A 0.28, and Doo ~ 0.25). As shown in Fig. 5 



the agreement with the numerical results is very good. 
Remarkably, we have obtained a similar Dq in the study 
of the quantum evolution matrix associated to certain 
classical interval-exchange maps whose level statistics is 
exactly given by SP 21] . Based on this finding we suggest 
that the Dq above describes the multifractal properties 
of all systems with SP. 

As a double check we have also evaluated Dq via the 
scaling of eigenfunction moments Vq = J d'^r\tp{r)\'^'^ cx 
j^-Dqiq-i) ^ff^iY^ respect to the sample size L = N. In par- 
ticular Dq was obtained by hnear regression of (logVq) 
with respect to logA^. As in the previous case, the av- 
erage is taken over all eigenvectors and different realiza- 
tions of the evolution and finite size effects are removed 
by an appropriate fitting. The Dq obtained in this way is 
in complete agreement with the previous one within the 
numerical errors. 

Finally we compare the above Dq with the predictions 
of critical statistics. First of all we remark that a compar- 
ison with a 3D (or 4D) disordered conductor at the AT 
is not entirely satisfactory since parameters as the slope 
of the number variance (^ 0.27 in the 3D case) defining 
the AT are different from the SP prediction (1/2). We 
thus find more appropriate to compare the SP results 
with those of a random banded matrix with a power-law 
(~ 1/r) decay and bandwidth b where the spectral 
statistics is given by critical statistics and the multifrac- 
tal character of the eigenstates has been established ana- 
lytically. In the limit 6^1, closer to SP, it is found that 
for q > 1/2 



D„ 



4b r(g-i/2) 



(7) 



As shown in Fig. 5 (inset), this function does not describe 
the multifractal dimensions associated to SP for any b. 
Thus we can conclude that despite their similarities criti- 
cal statistics and SP belong to different universality class. 



III. QUANTUM ANOMALOUS DIFFUSION 
AND EXPERIMENTAL VERIFICATION 



In this section we investigate how quantum transport 
properties are affected by the step-like non-analytical po- 
tential. Our motivation is twofold: On the one hand we 
wish to examine possible ways to test our results experi- 
mentally. On the other hand we are interested in finding 
out to what degree the phenomenon of dynamical local- 
ization typical of kicked rotors with smooth potentials 
^ cos{q) is affected by the classical singularity. We re- 
call that dynamical localization for smooth chaotic po- 
tentials manifests itself in the quantum suppression of 
classical diffusion in momentum space due to interfer- 
ence effects. Thus, contrary to the BGS conjecture [l|, 
though the classical dynamics is chaotic eigenstates are 
exponentially localized in momentum space and the level 
statistics is Poisson. 



8 




FIG. 6: (Color online) Classical and quantum variance {p'^{t)) 
of the density of probability P{p, t) (see text) for the Hamil- 
tonian Eq.lQ (T = 1,?!, = 1) and potential given by Eq.® 
with a — 10~*. The classical counterpart was obtained from 
10^ random initial conditions with zero initial momentum and 
uniformly distributed q in {~n,n). 
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FIG. 7: (Color online) Quantum coarse-grained P{k, t) versus 
k at various times t associated with the potential Eq. Q 
(«o = TT and a = n/2) plus a perturbation Vperiq) = 8cos(g). 
For a given pair of k and t, Pc(k, t) is the average of P{k' , t') 
over k' ^k- 63, fc -I- 64 and t' = t - 24, ...,t + 25. Pc{k,t) 
has power law tails ~ |fe|"^ for |fc| > t/hp where p 

is the spectral density . Similar results are observed if only 
the potential Eq. is considered. 



We start by computing both the quantum and classical 
density of probability P{k,t) = | (fc|0(t)) p of finding a 
particle with momentum p = kh at time t for a given 
initial state |0(O)) = |0) (initial condition in the classical 
case). 

The first problem we face is that the potential Eq.Q 
has a quite peculiar classical limit. The force exerted on 
the particle is a sum of two Dirac delta functions S{q±a), 
namely, it vanishes elsewhere except at q — ±a where 
formally it diverges. Consequently the classical motion 
is that of a free particle in a ring except for q = ±a 
where it gets an infinite force. Obviously in order to 
compare classical with quantum predictions one has first 
to smooth the step potential. Thus for the analysis of 



dynamical localization we consider the potential 

V{q) = Si{{a + q)/a) + Si{{a~q)/a) (8) 

where Si{q) — dt is the sine integral function. 

For CT ^ we recover the step-like potential. The 
classical force associated to this potential is F[q) — 
sin ^^—^ \- sin (Dirac delta functions in the 

a—q a a+q a v 

limit (7^0). With this potential we get a well de- 
fined classical limit for any finite a. We have computed 
the classical -P(p, t) by evolving the classical equation of 
motion for 10^ different random initial conditions with 
zero momentum p = and uniformly distributed po- 
sition along the interval (— tt, tt). We have found that 
P(p,t) ~ 21? for IpI < c{a)V2Dt where D - l/2cr 
and c((t) increases as a decreases. Outside this region 
P{p,t) resembles that of a standard diffusion process. 
Thus for sufficiently small p and t the diffusion is anoma- 
lous and then gradually becomes normal. In Fig. 6 we 
plot the quantum and classical variance {p^it)) of the 
density of probability P(p, t) {p = hk for the quantum 
case) measured at time t. Calculations were carried out 
for the differentiable potential of Eq.©. It is observed 
that, as expected, for short times both classical and quan- 
tum results coincide. However, after a certain breaking 
time tf, (X D 100 the quantum particle still diffuses 
but a slower rate than the classical one thus suggesting 
that diffusion is weaken by quantum interference effects. 
We relate this new region of weak dynamical localization 
the the effect of the classical singularity on the quantum 
dynamics. 

For even longer times ^ 20000 related to the 
crossover to normal classical diffusion the quantum par- 
ticle tends to localize in momentum space as a conse- 
quence of standard dynamical localization effects. The 
latter regime is due to the underlying smooth nature of 
the potential investigated. Thus in order to observe gen- 
uine quantum effects associated to classical singularities 
the value of a must be such that tt <C tc- We remark that 
this condition should be met by any experiment aiming 
to confirm the results reported in this letter. 

For the sake of completeness we have also investigated 
the specific form of the quantum P{k, t) as a function of 
k (from now on we switch back to our original potential 
Eq.®). In the context of a disordered conductor it has 
been reported that [23| at the AT the quantum difllusion 
is anomalous. For time scales large enough the density 
of probability has a power-law form with the exponent 
depending on the multifractal dimension D2. Similar 
results has also been recently obtained in momentum 
space for a kicked rotor with a logarithmic singularity 
[l6l |. As expected, we have also observed power-law tails 
in our model in both the short and the long t limit. For 
|fc| ^ t/Tip (with p = l/27r the spectral density) a best 
fit estimates yields P{k,t) ~ i0.84±o.oi|;j|-2 (ggg ^ig 7), 
similar to the classical prediction. By contrast in the op- 
posite limit though power-law tails have also been found 
the diffusion is slower in agreement with previous find- 
ings. 
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FIG. 8: (Color online) A typical eigenstate of the evolution 
matrix Eq.Q = 1, ^ = 1) with potential Eq. Q {vo — tv 
and a = n/2) plus a perturbation Vper ~ 8 cos q. The modulus 
of the eigenstate \tpk\ decays from its peak ko as a power-law 
with an exponent close to —1. Indeed the best fit (dashed 
line) corresponds to a slope —0.98. 



quantum localization. For the case of a smooth periodic 
standing wave, dynamical localization was also detected 
experimentally in full agreement with the theoretical pre- 
dictions. We thus propose that the full P{k, t), the accu- 
mulated probability or the variance represented in Fig. 6 
may by accessible to experimental verification provided 
that the smooth periodic standing wave is replaced by 
the step-like wave studied in this paper. Obviously from 
a experimental point of view the step-like singularity is 
only approximated. Typically the experimental signal is 
composed of a limited number of Fourier components and 
consequently it is smooth on sufficiently small scales (as 
the potential discussed in connexion with dynamical lo- 
calization). However, as we have shown in the previous 
section, we expect our results to hold in these 'almost' 
non-analytical potential at least up to a certain time scale 
related to the underlying smoothness of the potential. 



IV. CONCLUSIONS 



Another feature induced by the step-like singularity is 
the power-law localization of the eigenstates. In Fig. 8 
is clearly observed that the eigenstates have power- 
law tails with an exponent around minus one. This is in 
clear contrast with the exponential localization observed 
in the kicked rotor with a smooth potential. However is 
similar to the case of a potential with a log singularity 
[Tfil |: not surprisingly in both cases the off-diagonal ele- 
ments of the evolution operator present a similar power- 
law decay [l5l |. We mention there is no contradiction 
between this smooth power-law behavior and the multi- 
fractal features investigated previously. The point is that 
the multifractal character of the eigenstates appears as 
strong fluctuations around the smooth power-law behav- 
ior above. Roughly speaking we can say that smooth 
power-law localization is the precursor of multifractality. 

Finally we discuss how the findings of this paper may 
be tested experimentally. The first direct experimental 
realization of the quantum kicked rotor with a smooth 
potential was reported by Raizen |29| and coworkers in 
1995. The experimental set up consisted of a dilute 
sample of ultracold atoms (typically Cesium or Sodium) 
'kicked' by a periodic standing wave of near-resonant 
light that is pulsed on periodically in time to approxi- 
mate a series of delta functions. The typical output of 
the experiment is the distribution of the atom momen- 
tum as a function of time. It is also possible to estimate 
the quantum breaking time signaling the beginning of 



In this paper we have studied a kicked rotator with a 
non-analytical step-like singularity. It has been identified 
a region of parameters where the level statistics is exactly 
given by SP. These results are universal in the sense that 
they do not depend on the specific form of the potential 
but only on the presence of the classical singularity. The 
eigenfunctions have been shown to be multifractal but 
with a set of multifractal exponents Dq different from 
the prediction of critical statistics. We have also conjec- 
tured, based on the numerical analysis and the compari- 
son to other models with SP, that all systems described 
by SP have the same form of Dq {Dq = A/q + Doo)- Fi- 
nally we have studied transport properties. It has been 
found that, unlike the standard kicked rotor, dynamical 
localization slow down but does not stop quantum dif- 
fusion. We have also discussed the possibility of exper- 
imental verification of our findings by using ultra cold 
atoms kicked by a standing wave with an approximated 
step-like form. 
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